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We present a unified treatment of the prevalence of different double ionization (DI) pathways as a 
function of the sum of the final electron energies in strongly driven Helium. We do so as a function 
of laser frequency and intensity. At small total energy a new DI pathway prevails where both 
electrons are ionized with a delay after re-collision. We find that three-body collisions between the 
two electrons and the nucleus underly this pathway-we refer to it as triple collision (TC) pathway. 
At high total energy the pathway that prevails is the one where one electron is ionized with a delay 
following re-collision — we refer to it as Delayed pathway. Asymptotic observables that trace the 
prevalence of the latter pathways are identified. 

PACS numbers: 32.80.Rm, 31.90,+s, 32.80.Fb, 32.80.Wr 



I. INTRODUCTION 

The Helium atom driven by strong infrared laser fields 
has attracted considerable interest over the last two 
decades. The competing forces in strongly driven He, 
namely, the influence of the laser field on each electron, 
the electron-electron repulsion, and the interaction of 
each electron with the nucleus result in a wealth of phys- 
ical phenomena — many still unexplored. When the field 
intensity is large the two electrons are stripped out se- 
quentially [1 while for smaller intensities non-sequential 
double ionization (NSDI) dominates. In this work, we 
address the latter intensity regime. 

According to the three-step, also called re-scattering, 
model [2 — the accepted mechanism for NSDI — 1) one 
electron escapes through the field-lowered Coulomb po- 
tential, 2) it moves in the strong infrared laser field 
and 3) it returns to the core to transfer energy to the 
other electron remaining in He + . Theoretical support for 
the three step model has been provided, among others, 
by [3HS]. Although the re-scattering model has worked 
well in providing the interpretation of the basic strong 
field phenomena, such as ATI (Above Threshold Ion- 
ization), HHG (High-Order Harmonic Generation) and 
NSDI (Non-Sequential Double Ionization), recent kine- 
matically complete experimental investigations pj |H], 
have revealed additional structure in the latter for 800 
nm laser fields. Specifically they observed the so-called 
"finger-like" structure (V-shape) in the correlated mo- 
menta of the outgoing electrons, previously predicted by 
quantum mechanical studies of He DI at 400 nm [9 . Re- 
cent classical [lOl [TT] and quantum-mechanical [12] the- 
oretical treatments have added to our understanding of 
the latter. 

A still open problem is a unified description of the dif- 
ferent DI pathways the two electrons follow to escape 
after re-collision. The Direct and Delayed are two well 
established pathways: one involves an almost simulta- 
neous ionization (SI) of the two electrons and the other 



an ejection with a delay of approximately one quarter 
of a laser period or more (The delayed pathway is also 
referred to as re-collision-induced excitation with subse- 
quent field ionization, RESI [T3HT5] ). Taking a first step 
towards a general description of the DI pathways, in ref. 
[T6] we used full-dimensional classical and quantum me- 
chanical techniques for He at 400 nm. After establishing 
a remarkable agreement of the two treatments we used 
the classical calculation to explore the features of the DI 
mechanisms in three different regimes of the total elec- 
tron energy. We found that the delayed pathway prevails 
for small intensities — not addressed in the current work — 
independently of total electron energy. (Throughout this 
work the total electron energy we refer to is the sum of 
the final kinetic energies of the two electrons). In con- 
trast, at higher intensities the SI pathway prevails up to 
an upper-limit in total energy that shifts upwards with 
increasing intensity. The effect of the nucleus proved cru- 
cial for explaining this upward shift. 

In this work, in a quasiclassical framework, we explore 
the prevalence of the DI pathways as a function of total 
energy in energy steps as small as the immense compu- 
tational challenge of the endeavor allows — 600-700 thou- 
sand of DI events for the whole energy regime. We do 
so as a function of laser intensity and for two frequen- 
cies corresponding to wavelengths of 400 nm and 800 nm. 
We also identify asymptotic observables where the preva- 
lence of a DI pathway over the rest can be clearly traced. 
This general treatment on one hand allows us to explore 
whether previous findings at 400 nm are valid for 800 nm 
and on the other reveals new ones. 

One such new finding is that at very small total en- 
ergy a new DI pathway prevails. We refer to it as 
TC (triple collision) with both electrons ionizing more 
than a quarter of a laser cycle after re-collision. TC 
is present for all laser frequencies and intensities con- 
sidered in the current study. We find that three-body 
collisions between the two electrons and the nucleus un- 
derly the TC pathway. In addition, at high total en- 
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ergy we find that the Delayed pathway overtakes the Di- 
rect one during a very sharp transition. These findings 
hold true for a range of different intensities and frequen- 
cies, namely, for 9/13x 10 14 Watts/cm 2 at 400 nm and for 
2.25/3.25xl0 14 Watts/cm 2 at 800 nm. 

Besides the above mentioned features common to DI 
at 400 nm and 800 nm, we identify a frequency depen- 
dent feature. The number of times the re-colliding elec- 
tron returns to the core before both electrons are ion- 
ized decreases with increasing frequency. (This num- 
ber is defined through electron returns roughly at times 
(2/3+n)T — three-step model, with T the period of the 
field). This difference roughly accounts for the differ- 
ent patterns of the correlated momenta for the Delayed 
mechanism observed at 400 nm and 800 nm. 



II. MODEL 

The three-dimensional quasiclassical model we use en- 
tails one electron tunneling through the field- lowered- 
Coulomb potential with a quantum tunneling rate given 
by the ADK formula [T?] . The laser pulse used in the 
classical calculations is E(t) = Eo(t) cos(ut) and is lin- 
early polarized along the z-axis. The longitudinal mo- 
mentum is zero while the transverse one is given by a 
Gaussian distribution. The remaining electron is mod- 
eled by a microcanonical distribution [18]. The pulse 
envelope is defined as Eo(t) = Eq (a constant) for 
< t < 6T and E (t) = E cos 2 (u(t - 6T)/12) for 
6T < t < 9T with T the period of the field. An ad- 
vantage of our classical propagation over other classical 
techniques is that we employ regularized coordinates [19] 
(to account for the Coulomb singularity) which results in 
a faster and more stable numerical propagation. 



III. DOUBLE IONIZATION MECHANISMS 

A. % Contribution of DI mechanisms as a function 
of total electron energy 

In Fig. [T] we show the percentage contribution of 
the DI mechanisms as a function of the sum of the 
two electron energies. We do so for laser pulse inten- 
sities of 2.25/3.25 xl0 14 Watts/cm 2 at 800 nm and for 
9/13xl0 14 Watts/cm 2 at 400 nm. The range of laser in- 
tensities we currently consider is important because, at 
400 nm, for 9 (12) xlO 14 W/cm 2 the maximum return 
energy of the re-colliding electron (3.2 U p according to 
the three-step model [2]) equals the first excitation (ion- 
ization) energy of ground state He + and similarly for 2.25 
(3.25) xlO 14 W/cm 2 at 800 nm. The ponderomotive en- 
ergy Up = Eq/(4uj 2 ) is the average energy an electron 
gains in one laser cycle. The intensities 9 (13) xlO 14 
W/cm 2 at 400 nm correspond to the same ponderomo- 
tive energy as for 2.25 (3.25) xlO 14 W/cm 2 at 800 nm, 
and thus the same Keldysh parameter 7 = \/I p / (2U P ) 



[20] . We adapt an energy step, for the final total energy, 
of 0.4 Up to achieve good statistics in all energy intervals. 
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FIG. 1. First row: probability distribution for 61+62; Middle 
row: % of DI mechanisms with • for SI, ♦ for RESIa, * for 
RESIb and A for TC; Bottom row: energy sharing, i.e., r~^~j 
as a function of ei + 62- (From left to right) first column: 
2.25 xlO 14 Watts/cm 2 at 800 nm, second column: 3.25 xlO 14 
Watts/cm 2 at 800 nm, third column: 9xl0 14 Watts/cm 2 at 
400 nm and fourth column: 13xl0 14 Watts/cm 2 at 400 nm. 

To identify the main DI energy transfer pathways we 
use the time delay between the re-collision of the free 
electron with the parent ion and the onset of ionization 
of the second electron [21]. For a definition of the time 
of ionization see [10] and references there in. In the di- 
rect ionization pathway (SI) both electrons are ionized 
simultaneously very close (less than a quarter of a laser 
period) to the re-collision time. In the delayed ionization 
pathway, the re-colliding electron excites the remaining 
electron but does not ionize it. The electron is subse- 
quently ionized at a peak (RESIa) or at a zero (RESIb) 
of the laser electric field [I3j [14] . Finally in the TC path- 
way both electrons are ionized more than a quarter of a 
laser period following re-collision. The % break-up of the 
different DI mechanisms as a function of the total elec- 
tron energy exposes universal features present for differ- 
ent laser intensities and frequencies. 

B. A new pathway of DI ionization for small total 
electron energy 

One such feature is that for small total energy the two 
electrons escape through the TC pathway; both electrons 
are ionized with a delay of more than a quarter of a laser 
cycle with respect to the re-collision time. TC reaches 
already 86% for 2.25xl0 14 Watts/cm 2 at 800 nm in the 
energy interval [0, 0.4] U p , making this parameter regime 
ideal for the study that follows. 
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The signature feature of the TC pathway is three-body 
collisions between the nucleus and the two electrons. For 
the vast majority of trajectories the re-colliding electron 
following re-collision gets temporarily "trapped" orbiting 
the nucleus. In Fig. [2] a) a representative TC trajectory 
is depicted where the first electron, following re-collision, 
orbits the nucleus once before both electrons are ionized. 
Trajectories where the first electron orbits the nucleus 
more than once are also common. The three body colli- 
sion is best evidenced by the subset of trajectories where, 
in addition to orbiting the nucleus, the two electrons ion- 
ize at the same time, see Fig. [2] if the interaction of the 
laser field with the two electrons was the prevalent one 
then the electronic positions along the field polarization 
would be similar. Instead, the three-body collision pre- 
vails and since the energy shared among the three parti- 
cles is small the two electrons escape in almost opposite 
directions along the field polarization. 

To better understand why the first electron gets 
trapped after re-collision takes place, we identify the 
phase of the laser field when the first electron tunnels. 
In Fig. U we show that in the TC pathway the re- 
colliding electron tunnels into the continuum when the 
phase of the laser field is around -3°. This results in 
the re-colliding electron returning to the nucleus with 
small energy compared to the maximum return energy 
of 3.2 Up achieved when tunneling takes place around 
17° (three-step model). In contrast, in the SI mechanism 
the re-colliding electron tunnels when the phase of the 
field is around 20°, resulting in higher energy at the time 
of re-collision. The return to the nucleus with small en- 
ergy seems to be a necessary condition for the signature 
feature of the TC pathway. We finally note that for the 
majority of the TC trajectories re-collision takes place 
half a cycle before the re-colliding electron ionizes and 
the second electron is ionized primarily at the same time 
as the first electron or half a cycle later. 
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FIG. 2. a) The radii of the re-colliding (black line) and of 
the initially bound electron (grey line); b) The distances in 
the direction of the laser polarization (z) for the re-colliding 
(black dashed line) and for the initially bound electron (grey 
dashed line) for TC for 2.25xl0 14 Watts/cm 2 at 800 nm. The 
re-collision time is around 1.75 laser cycles. Both electrons are 
ionized between 2-2.5 cycles. 

The strong electron-electron repulsion is further evi- 
denced using a triple differential distribution in Fig. [4] 
For the smallest energy interval considered we compute 
the DI probability as a function of the sum of the mo- 




phase of the laser field 



FIG. 3. Probability distribution for the re-colliding electron 
to tunnel at a given phase of the laser field; • for SI, ♦ for 
RESIa, ★ for RESIb and < for TC for 800 nm and 2.25xl0 14 
Watts/cm 2 . 
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Pl,z+P2,z( a - U -) 

FIG. 4. DI probability as a function of the sum of the x 
component of the momenta versus the z component of the 
momenta, for 2.25xl0 14 Watts/cm 2 at 800 nm for the four 
different mechanisms SI (left), RESIa (second from left) RE- 
SIb (third from left) TC (right) for ei + e 2 in the [0,0.4] U p 
energy interval. 



menta components parallel to the field polarization (z- 
axis) and perpendicular to it. From Fig.JIJd) we find that 
the highest DI probability for the TC pathway is when 
both the sum of the parallel and perpendicular compo- 
nents of the momenta are zero suggesting that a strong 
electron-electron repulsion underlies the two electron es- 
cape. This is not the case for the other pathways as 
shown in Fig. [4] a), b), c). Thus, one way to confirm ex- 
perimentally the prevalence of the TC pathway for small 
energies would be to measure the triple differential prob- 
ability described above in different energy regimes and 
find that only for the smaller energy interval the triple 
differential distribution resembles the TC one in Fig. [4] 
d). Indeed, in Fig. [5] we plot the correlated sum of the 
momenta in four energy regimes accounting for all dou- 
bly ionizing trajectories (all pathways included) and find 
that only for the lowest energy interval, see Fig. [5] a), the 
correlated sum of the momenta resembles that of Fig. [4] 

d). 

We have provided evidence that each electron is 
trapped by the nucleus after re-collision (Fig. thus 
a strong electron-nucleus interaction takes place and 
we have shown that after re-collision there is a strong 
electron-electron repulsion (Fig. [4j Fig.|5|. We thus con- 
jecture that triple collisions — interaction of the two elec- 
trons with the nucleus — underlies the pathway prevail- 
ing for small total electron energy. The TC pathway is a 
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FIG. 5. DI probability as a function of the sum of the x 
component of the momenta versus the z component of the 
momenta, for 2.25xl0 14 Watts/cm 2 at 800 nm for four en- 
ergy intervals; from left to right [0,0.4] U p , [0.8,1.2] U p , [3.6, 
4.0] Up, [5.6, 6.0] Up. The sum of the momenta is scaled by 
l/\/Emid, where E mi d is the middle of the respective energy 
interval. 



trace of collisional physics in strongly driven He, since a 
three-body interaction between two electrons and a nu- 
cleus dominates the threshold behavior of electron impact 
and photo-ionization processes, see Wannier's law [23] . 

Finally, let us note that the three-body collisions are 
more prevalent the smaller the interaction with the field 
is, explaining the decrease in the % contribution of the 
TC pathway with increasing frequency and intensity, see 
Fig. [I] second row. 



C. Direct versus Delayed pathways 

From Fig. [5] we see that as the total energy increases 
the correlated sum of the momenta changes significantly 
indicating the prevalence of DI pathways other than the 
TC one. This is the focus of the current section. 

From Fig. [T] we see that as the energy increases the 
Delayed pathway prevails up to total energy that differs 
with laser frequency and intensity; for larger intensities it 
extends up to around 1U P while for 2.25 x 10 14 Watts/cm 2 
at 800 nm up to 3U P . The SI mechanism prevails for in- 
termediate total energies extending up to approximately 
5 Up for small intensities and up to around 7-8 U p for 
higher ones, see Fig.[l]and [16 j. This increase with inten- 
sity in the SI upper energy limit was attributed to more 
effective encounters of the re-colliding electron with the 
nucleus [16]. The prevalence of the SI pathway is more 
pronounced for higher intensities reaching almost 70%, 
see Fig. [I] Increasing further the total energy we find 
that, through a sharp transition, the Delayed pathway 
prevails over the Direct one. The highest total energy we 
consider is 8 XJ P in order to ensure good statistics in all 
energy intervals. 

The effect of the nucleus and the time of re-collision 
can be extracted from the classical trajectories by plot- 
ting for each DI pathway the average value of each elec- 
tron's potential energy — 2/|ri/ 2 | and momenta compo- 
nents parallel to the field polarization. As for 400 nm 
[16], we find that the re-collision time at 800 nm shifts 
from 0.5 T (T is the laser period) at small total energies 
to (n+2/3)T (three-step model), n = 0, 1, 2, 3, ... at high 
total energies. The effect of the nucleus at 800 nm as 



for 400 nm is very important for total energies above ap- 
proximately 5 Up. However, the effect of the nucleus is 
also crucial for very small total energy as well where the 
TC pathway prevails, as we have shown above. 

Another trace of collisional physics on strong field ion- 
ization is the overall shape of the SI contribution to DI 
as a function of total energy. The latter is very small for 
small and high energies while it is large for intermediate 
ones. Thus, SI qualitatively resembles electron impact 
ionization of He + as a function of excess energy (energy 
above the ionization threshold of He + ) [22]. While in the 
field-free case the impacting electron can have any initial 
energy, for the driven system the maximum return en- 
ergy to the core is 3.2 \J P (three-step model), defining a 
different energy scale for the two processes. 

Given the above, for high total final energy, in the SI 
pathway, the impacting electron does not transfer enough 
energy at the time of re-collision to the second electron 
to ionize both electrons. However, the presence of the 
field makes possible a transfer of energy to the second 
electron even after re-collision through the RESI pathway 
(Delayed pathway) thus casting the RESI pathway the 
prevailing one for high total energy. 

Asymptotic Observables. The double differential prob- 
ability of the energy sharing |ei — 62 1/ (ei + e 2 ) versus the 
total electron energy reveals additional information for 
the DI process. For all four intensities currently consid- 
ered we find that for increasing total energy the two elec- 
trons share the energy more asymmetrically compared 
to smaller energies. An asymmetry in energy sharing is 
also observed in electron impact ionization for high ex- 
cess energy — many times referred to as U-shape energy 
sharing [24] . When the Delayed pathway prevails at high 
total energies, the energy sharing peaks around approx- 
imately 0.4-0.6, see Fig. [I] third row. A study of the 
double energy differentials in energy (not shown here) 
reveals that the initially bound electron escapes with an 
energy around 2 U p (the energy a "free" electron gains 
in the laser field). The re-colliding electron escapes with 
almost 3 times larger energy than the initially bound 
electron — the asymmetry is larger for higher intensities 
due to stronger backscattering of the first electron re- 
sulting in its faster escape. Thus, this double differential 
probability can serve as a good predictor of the value 
of the total electron energy (for high energy) where the 
Delayed mechanism prevails. 

RES la versus RE Sib. In Table I, we show a break-up 
of the % contribution of the DI mechanisms independent 
of energy. We also show results for 8.12 x 10 13 Watts/cm 2 
at 1600 nm corresponding to the same Keldysh parame- 
ter as that for 3.25/13 xlO 14 Watts/cm 2 at 800/400 nm. 
Table I shows that when the Keldysh parameter is kept 
constant, smaller frequencies favor the RESIa versus the 
RESIb meaning that the initially bound electron ionizes 
mostly around a maximum and not a zero of the field. 
This is to be expected since the larger the frequency the 
faster the velocity changes around a zero of the field re- 
sulting in higher contribution of the RESIb pathway. 
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TABLE I. % of DI mechanisms over all DI events. 

For each intensity, all mechanisms listed sum up to « 100% 



8.12 xlO 13 W/cm 2 at 1600 nm 


SI RESIa RESIb DE 




54.7 27.2 16.5 0.8 


2.25 xlO 14 W/cm 2 at 800 nm 


SI RESIa RESIb DE 




26.7 24.95 40.19 7 


3.25 xlO 14 W/cm 2 at 800 nm 


SI RESIa RESIb DE 




60.8 18.6 18.1 1.8 


9 xlO 14 W/cm 2 at 400 nm 


SI RESIa RESIb DE 




45.2 13.1 35.3 5.5 


13 xlO 14 W/cm 2 at 400 nm 


SI RESIa RESIb DE 




53.6 15.41 26.5 3.53 



3/8xl0- 4 for 2.25/3.25xl0 14 Watts/cm 2 at 800 nm, and 
2/3 xlO" 3 for 9/13xl0 14 Watts/cm 2 at 400 nm. 

We find that this difference in the number of re- 
collisions affects the correlated momenta of the RESI 
pathway. Specifically, we see from Fig. [7] that DI through 
the RESIa pathway at 400 nm favors cross-shaped cor- 
related momenta along the two axes with the two elec- 
trons sharing the energy asymmetrically. At 800nm, DI 
through the RESIa pathway favors a structure in all four 
quadrants — as for 400 nm — however, the two electrons 
share the energy more symmetrically. A break down of 
the RESIa correlated momenta to contributions from dif- 
ferent number of returns to the core of the re-colliding 
electron reveals that DI after one re-collision favors the 
cross-shaped structure along the two axes. DI following 
more than one re-collisions favors a superposition of the 
two structures discussed above. Thus the prevalence of 
different number of re-collisions, before DI takes place, 
for different frequencies roughly explains why our RESIa 
correlated momenta are different for 400 nm and 800 nm. 
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D. Number of re-collisions and transfer of energy 

We have found that for small total energy — where the 
TC pathway prevails — the transfer of energy from the re- 
colliding to the bound electron takes place first through 
a re-collision and subsequently through a triple collision. 
For the Direct and the Delayed pathway which have over- 
all (independent of total energy) the largest contribu- 
tion to DI the transfer of energy takes place through re- 
collisions. 

A legitimate question is how many re-collisions are re- 
quired before the first electron transfers sufficient en- 
ergy to the second one so that they both ionize and 
whether this number changes with frequency. Fig. [6] 
shows the number of re-collisions required for DI for the 
RESI pathway but similar results hold for the SI path- 
way as well. We find that at 400 nm the two elec- 
trons doubly ionize mostly after one re-collision with 
the second re-collision also contributing significantly to 
DI. The second re-collision has the largest contribution 
to DI at 800 nm. Further decreasing the frequency 
to 1600 nm, we find that the first re-collision's contri- 
bution to DI is only minor while the contribution of 
the third one is as prominent as the second's. The 
decreasing number of re-collisions with increasing fre- 
quency suggests a more efficient transfer of energy to 
the initially bound electron for larger frequencies. In- 
deed, the DI probability increases with increasing fre- 
quency: 2.4 xlO" 4 for 8.12xl0 13 Watts/cm 2 at 1600 nm, 



FIG. 6. The % contribution of different number of re- 
collisions for the RESI mechanism out of all DI events. • 
for 1 re-collision, ★ for 2 re-collisions, ■ for 3 re-collisions; 
Point A/B corresponds to 9/13x 10 14 Watts/cm 2 at 400 nm, 
C/D to 2.25/3.25xl0 14 Watts/cm 2 at 800 nm, and E to 
8.12xl0 13 Watts/cm 2 at 1600 nm. 




pja.u.) 

FIG. 7. The correlated momenta for the RESIa mechanism 
(a) and RESIb mechanism (b); 1) 2.25 xlO 14 Watts/cm 2 at 
800 nm , 2) 9xl0 14 Watts/cm 2 at 400 nm. 



IV. CONCLUSIONS 

We have shown that the prevalence of the Direct and 
Delayed pathway for intermediate and high total energy, 
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respectively, is common for laser pulses at 400 nm and the correlated sum of the momenta components. 
800 nm for a range of intensities. We have also identified 
TC as a new DI pathway that prevails for small total en- 
ergy. Three-body collisions between the nucleus and the 
two electrons underly the TC pathway. Given the small 
energy available to the system the two electrons escape 
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